Plot and fit distributions of variables

Histogram of Velocity of all

Log distribution, 0’s make up 3% of all data

CC.TotalData <- na.omit(TotalData)
head(CC.TotalData)
10^mean(log10(CC.TotalData$v[CC.TotalData$v>0]))
[1] 0.0348809
sd(log10(CC.TotalData$v[CC.TotalData$v>0]))
[1] 0.5467087

Now we’re going to plot each unique krill and their bimodality with a dip test (less than 0.05 is multimodality)

Starting to look at bimodality in the variables by merging the TotalData frame with the “tab” table

Now looking at turning angles

for (i in 1:length(ind)){
plot(CC.TotalData$turn.angle[CC.TotalData$D_V_T==ind[i]], log10(CC.TotalData$v[CC.TotalData$D_V_T==ind[i]]),
            xlab = "Turn Angles",
     ylab = "Velocity",
          main = "") 
}

NA

save.image("~/Post-doc/Data/Total Merged Data File.RData")
LS0tDQp0aXRsZTogIkNoZWNraW5nIERpc3RyaWJ1dGlvbiBmb3IgVmFyaWFibGVzIg0Kb3V0cHV0OiBodG1sX25vdGVib29rDQotLS0NCg0KUGxvdCBhbmQgZml0IGRpc3RyaWJ1dGlvbnMgb2YgdmFyaWFibGVzDQoNCmBgYHtyfQ0KDQpsb2FkKCJDOlxcVXNlcnNcXG5pY29sZWgzXFxEb2N1bWVudHNcXFBvc3QtZG9jXFxEYXRhXFxUb3RhbCBEYXRhIE1lcmdlZCBGaWxlIChTZXAgMjcgMjAyMSkuUmRhIikNCg0KaGlzdChsb2cxMChUb3RhbERhdGEkdiksDQogICAgIHhsYWIgPSAiVmVsb2NpdHkgKExvZ14xMCBtbS9zKSIpDQoNCmBgYA0KDQpIaXN0b2dyYW0gb2YgVmVsb2NpdHkgb2YgYWxsDQoNCkxvZyBkaXN0cmlidXRpb24sIDAncyBtYWtlIHVwIDMlIG9mIGFsbCBkYXRhDQoNCmBgYHtyfQ0KQ0MuVG90YWxEYXRhIDwtIG5hLm9taXQoVG90YWxEYXRhKQ0KaGVhZChDQy5Ub3RhbERhdGEpDQoxMF5tZWFuKGxvZzEwKENDLlRvdGFsRGF0YSR2W0NDLlRvdGFsRGF0YSR2PjBdKSkNCnNkKGxvZzEwKENDLlRvdGFsRGF0YSR2W0NDLlRvdGFsRGF0YSR2PjBdKSkNCg0KYGBgDQpOb3cgd2UncmUgZ29pbmcgdG8gcGxvdCBlYWNoIHVuaXF1ZSBrcmlsbCBhbmQgdGhlaXIgYmltb2RhbGl0eSB3aXRoIGEgZGlwIHRlc3QgKGxlc3MgdGhhbiAwLjA1IGlzIG11bHRpbW9kYWxpdHkpDQoNCmBgYHtyfQ0KaW5kIDwtIHVuaXF1ZShUb3RhbERhdGEkRF9WX1QpDQpwcmludChpbmQpDQpsZW5ndGgoaW5kKQ0KbGlicmFyeShkaXB0ZXN0KQ0KbGlicmFyeShEZXNjVG9vbHMpDQoNCnRhYiA8LSBtYXRyaXgoZGF0YSA9IE5BLCBucm93ID0gMTM1LCBuY29sID01LCBieXJvdyA9IFQpDQpjb2xuYW1lcyh0YWIpIDwtIGMoJ2RpcC50ZXN0JywgJ3NrZXcnLCAnbWVhbi52ZWxvY2l0eScsICdzZC52ZWxvY2l0eScsICdJbmQnKQ0KDQp0YWIgPC0gYXMuZGF0YS5mcmFtZSh0YWIpDQp0YWIkSW5kIDwtIGluZA0KbWVhbi52ZWwgPC0gTlVMTA0Kc2QudmVsIDwtIE5VTEwNCmQudiA8LSBOVUxMDQpzLnYgPC0gTlVMTA0KDQpmb3IgKGkgaW4gMTpsZW5ndGgoaW5kKSl7DQptZWFuLnYgPC0gbWVhbihDQy5Ub3RhbERhdGEkdltDQy5Ub3RhbERhdGEkRF9WX1Q9PWluZFtpXV0pDQpzZC52IDwtIHNkKENDLlRvdGFsRGF0YSR2W0NDLlRvdGFsRGF0YSREX1ZfVD09aW5kW2ldXSkNCg0KbWVhbi52ZWwgPC0gcmJpbmQobWVhbi52ZWwsIG1lYW4udikNCnNkLnZlbCA8LSByYmluZChzZC52ZWwsIHNkLnYpDQoNCg0KdmVscyA8LSAoQ0MuVG90YWxEYXRhJHZbQ0MuVG90YWxEYXRhJERfVl9UPT1pbmRbaV1dKQ0KdmVscyA8LSBsb2cxMCh2ZWxzW3ZlbHM+MF0pDQpkIDwtIGRpcC50ZXN0KHZlbHMpDQpkLnAgPC0gZCRwLnZhbHVlDQpkLnYgPC0gcmJpbmQoZC52LCBkLnApDQoNCnMgPC0gU2tldyh2ZWxzKQ0Kcy52IDwtIHJiaW5kKHMudiwgcykNCiMjfQ0KDQojI2luZCA8LSB1bmlxdWUoQ0MuVG90YWxEYXRhJERfVl9UKQ0KDQojI2ZvciAoaSBpbiAxOmxlbmd0aChpbmQpKXsNCiBoaXN0KGxvZzEwKENDLlRvdGFsRGF0YSR2W0NDLlRvdGFsRGF0YSREX1ZfVD09aW5kW2ldXSksDQogICAgIGJyZWFrcyA9IDUwLA0KICAgICB4bGFiID0gIlZlbG9jaXR5IChMb2deMTAgbW0vcykiLA0KICAgICBtYWluID0gaW5kW2ldLA0KICAgICBzdWIgPSBkLnApICAjIyBjaGFuZ2UgdG8gZC5wIG9yIHMgdG8gcHJpbnQgdGhlIGRpcCB0ZXN0IG9yIHNrZXcgdmFsdWUgYXMgdGhlIHRpdGxlIGluc3RlYWQNCn0NCg0KdGFiJHNrZXcgPC0gcy52DQp0YWIkZGlwLnRlc3QgPC0gZC52DQp0YWIkbWVhbi52ZWxvY2l0eSA8LSAobG9nMTAobWVhbi52ZWwpKQ0KdGFiJHNkLnZlbG9jaXR5IDwtIChsb2cxMChzZC52ZWwpKQ0KdGFiDQp3cml0ZS50YWJsZSh0YWIsIGZpbGUgPSAifi9Qb3N0LWRvYy9EYXRhL2RpcC50ZXN0LnNrZXcudmVscy5jc3YiLCBzZXAgPSAiLCIsIGNvbC5uYW1lcyA9IFRSVUUpDQoNCnBsb3QodGFiKQ0KDQoNCmhlYWQoVG90YWxEYXRhKQ0Kc3RyKFRvdGFsRGF0YSkNClRvdGFsRGF0YSRGbG93LnJhdGUgPC0gYXMuZmFjdG9yKFRvdGFsRGF0YSRGbG93LnJhdGUpDQpUb3RhbERhdGEkQ2hsb3JvcGh5bGwgPC0gYXMuZmFjdG9yKFRvdGFsRGF0YSRDaGxvcm9waHlsbCkNCg0KbGlicmFyeShnZ3Bsb3QyKQ0KZ2dwbG90KFRvdGFsRGF0YSxhZXMoeD1GbG93LnJhdGUsIHk9dHVybi5hbmdsZSwgZmlsbD1DaGxvcm9waHlsbCkpKw0KICBnZW9tX2JveHBsb3Qobm90Y2g9Riwgbm90Y2h3aWR0aD0wLjMsb3V0bGllci5zaGFwZT0xLG91dGxpZXIuc2l6ZT0yLCBjb2VmPTEuNSkrDQogIHRoZW1lKGF4aXMudGV4dD1lbGVtZW50X3RleHQoY29sb3I9ImJsYWNrIikpKw0KICB0aGVtZShheGlzLnRleHQueD1lbGVtZW50X3RleHQoYW5nbGU9OTAsaGp1c3Q9MSx2anVzdD0wLjQpKSsNCiAgdGhlbWUocGFuZWwuZ3JpZC5taW5vcj1lbGVtZW50X2JsYW5rKCkpKw0KICBsYWJzKHNpemU9ICIiLHggPSAiRmxvdyBSYXRlIChjbS9zKSIsIHkgPSAiVHVybiBhbmdsZSAoZGVncmVlcykiLCB0aXRsZSA9ICIgICAgICAgICAgICAgICBMaWdodCIpKw0KICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXM9YygiZ3JlZW55ZWxsb3ciLCAiZ3JlZW4iLCJncmVlbjEiLCAiZ3JlZW4yIiwgImdyZWVuMyIsICJncmVlbjQiKSxuYW1lID0gIkNobG9yb3BoeWxsIChtZy9MKSIsDQogICAgICAgICAgICAgICAgICAgIGxhYmVscz1jKCIwIiwgIjQuNiIsICI2LjEiLCAiNy42IiwgIjEzLjUiLCAiMTkiKSkrDQogIGZhY2V0X2dyaWQofkxpZ2h0LCBzY2FsZXMgPSAiZnJlZV94Iiwgc3BhY2UgPSAiZnJlZSIpDQpgYGANCg0KDQpTdGFydGluZyB0byBsb29rIGF0IGJpbW9kYWxpdHkgaW4gdGhlIHZhcmlhYmxlcyBieSBtZXJnaW5nIHRoZSBUb3RhbERhdGEgZnJhbWUgd2l0aCB0aGUgInRhYiIgdGFibGUNCmBgYHtyfQ0KIyNhZ3JlZ2F0aW5nIENvbXBsZXRlIGNhc2VzIG9mIFRvdGFsRGF0YSBzbyB3ZSBjYW4gbWVyZ2UgaXQgd2l0aCB0YWIgZGF0YQ0KDQpBR0dfVEQgPC0gYWdncmVnYXRlKFRvdGFsRGF0YSwgYnkgPSBsaXN0KFRvdGFsRGF0YSREX1ZfVCksIEZVTiA9IG1lYW4pDQpoZWFkKEFHR19URCkNCkFHR19URCA8LSBBR0dfVERbIC1jKDI6MywgODoxMCkgXQ0KY29sbmFtZXMoQUdHX1REKSA8LSBjKCJJbmQiLCAiWCIsICJZIiwgIloiLCAiVHJhY2siLCAiRmxvdy5SYXRlIiwgIkNobG9yb3BoeWxsIiwgIkd1YW5vIiwgIkxpZ2h0IiwgImR4IiwgImR5IiwgImR6IiwgImQiLCAidngiLCAidnkiLCAidnoiLCAidiIsICJoZWFkaW5nIiwgInBpdGNoIiwgInR1cm5hbmdsZS54eSIsICJ0dXJuYW5nbGUueXoiLCAidHVybi5hbmdsZSIpDQpoZWFkKEFHR19URCkNCnRhaWwoQUdHX1REKQ0KDQoNCnRhYl9BR0cgPC0gbWVyZ2UoQUdHX1RELCB0YWIsIGJ5ID0gIkluZCIpDQpoZWFkKHRhYl9BR0cpDQpzdHIodGFiX0FHRykNCnRhYl9BR0ckRmxvdy5SYXRlIDwtIGFzLmZhY3Rvcih0YWJfQUdHJEZsb3cuUmF0ZSkNCnRhYl9BR0ckQ2hsb3JvcGh5bGwgPC0gYXMuZmFjdG9yKHRhYl9BR0ckQ2hsb3JvcGh5bGwpDQp0YWJfQUdHJEd1YW5vIDwtIGFzLmZhY3Rvcih0YWJfQUdHJEd1YW5vKQ0KdGFiX0FHRyRMaWdodCA8LSBhcy5mYWN0b3IodGFiX0FHRyRMaWdodCkNCg0KDQpmcmVxIDwtIHRhYmxlKHRhYl9BR0ckRmxvdy5SYXRlLCB0YWJfQUdHJENobG9yb3BoeWxsKQ0KcHJpbnQoZnJlcSkNCnByb2IgPC0gcHJvcC50YWJsZShmcmVxKSAjI1JlbGF0aXZlIEZyZXF1ZW5jeSBUYWJsZQ0KcHJpbnQgKHByb2IpDQoNCiAjI3N0YXJ0aW5nIHRvIHBsb3QgdGhlIGRpcCB0ZXN0IGFuZCBza2V3IGluIHRoZSB2YXJpYWJsZXMNCnBsb3QodGFiX0FHRyRGbG93LlJhdGUsIHRhYl9BR0ckZGlwLnRlc3QsIHhsYWIgPSAiRmxvdyBSYXRlIChjbS9zKSIsIHlsYWIgPSAiRGlwIFRlc3QgKHAudmFsdWUpIikNCnBsb3QodGFiX0FHRyRGbG93LlJhdGUsIHRhYl9BR0ckc2tldywgeGxhYiA9ICJGbG93IFJhdGUgKGNtL3MpIiwgeWxhYiA9ICJTa2V3IFRlc3QgKHAudmFsdWUpIikNCnBsb3QodGFiX0FHRyRDaGxvcm9waHlsbCwgdGFiX0FHRyRkaXAudGVzdCwgeGxhYiA9ICJDaGxvcm9waHlsbCAobWcvTCkiLCB5bGFiID0gIkRpcCBUZXN0IChwLnZhbHVlKSIpDQpwbG90KHRhYl9BR0ckQ2hsb3JvcGh5bGwsIHRhYl9BR0ckc2tldywgeGxhYiA9ICJDaGxvcm9waHlsbCAobWcvTCkiLCB5bGFiID0gIlNrZXcgVGVzdCAocC52YWx1ZSkiKQ0KDQpsaWJyYXJ5KGdncGxvdDIpDQpnZ3Bsb3QodGFiX0FHRyxhZXMoeD1GbG93LlJhdGUsIHk9c2tldywgZmlsbD1DaGxvcm9waHlsbCkpKw0KICBnZW9tX2JveHBsb3Qobm90Y2g9Riwgbm90Y2h3aWR0aD0wLjUsb3V0bGllci5zaGFwZT0xLG91dGxpZXIuc2l6ZT0yLCBjb2VmPTEuNSkrDQogIHRoZW1lKGF4aXMudGV4dD1lbGVtZW50X3RleHQoY29sb3I9ImJsYWNrIikpKw0KICB0aGVtZShheGlzLnRleHQueD1lbGVtZW50X3RleHQoYW5nbGU9OTAsaGp1c3Q9MSx2anVzdD0wLjQpKSsNCiAgdGhlbWUocGFuZWwuZ3JpZC5taW5vcj1lbGVtZW50X2JsYW5rKCkpKw0KICBsYWJzKHNpemU9ICIiLHggPSAiRmxvdyBSYXRlIChjbS9zKSIsIHkgPSAiU2tldyAocCB2YWx1ZSkiLCB0aXRsZSA9ICIiKSsNCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzPWMoImdyZWVueWVsbG93IiwgImdyZWVuIiwiZ3JlZW4xIiwgImdyZWVuMiIsICJncmVlbjMiLCAiZ3JlZW40IiksbmFtZSA9ICJDaGxvcm9waHlsbCAobWcvTCkiLA0KICAgICAgICAgICAgICAgICAgICBsYWJlbHM9YygiMCIsICI0LjYiLCAiNi4xIiwgIjcuNiIsICIxMy41IiwgIjE5IikpICANCiAjIyArIGZhY2V0X2dyaWQofkxpZ2h0LCBzY2FsZXMgPSAiZnJlZV94Iiwgc3BhY2UgPSAiZnJlZSIpDQoNCmBgYA0KDQpOb3cgbG9va2luZyBhdCB0dXJuaW5nIGFuZ2xlcw0KDQpgYGB7cn0NClRvdGFsRGF0YSR0dXJuLmFuZ2xleHkgPC0gYXRhbjIoVG90YWxEYXRhJFgsIFRvdGFsRGF0YSRZKQ0KVG90YWxEYXRhJHR1cm4uYW5nbGV5eiA8LSBhdGFuMihUb3RhbERhdGEkWSwgVG90YWxEYXRhJFopDQoNCg0KbHRoIDwtIGRpbShUb3RhbERhdGEpWzFdDQpkeDEgPC0gVG90YWxEYXRhJGR4WzE6KGx0aC0xKV0NCmR4MiA8LSBUb3RhbERhdGEkZHhbMjpsdGhdDQpkeTEgPC0gVG90YWxEYXRhJGR5WzE6KGx0aC0xKV0NCmR5MiA8LSBUb3RhbERhdGEkZHlbMjpsdGhdDQpkejEgPC0gVG90YWxEYXRhJGR6WzE6KGx0aC0xKV0NCmR6MiA8LSBUb3RhbERhdGEkZHpbMjpsdGhdDQpEIDwtIChkeDEqZHgyKSsoZHkxKmR5MikrKGR6MSpkejIpDQpkMSA8LSBzcXJ0KGR4MV4yICsgZHkxXjIgK2R6MV4yKQ0KZDIgPC0gc3FydChkeDJeMiArIGR5Ml4yICtkejJeMikNCg0KZGQgPC0gRC9kMS9kMg0KaGlzdChhY29zKGRkKS9waSoxODApDQoNClRvdGFsRGF0YSR0dXJuLmFuZ2xlIDwtIGMoTkEsIGFjb3MoRC9kMS9kMikpL3BpKjE4MA0KQ0MuVG90YWxEYXRhIDwtIG5hLm9taXQoVG90YWxEYXRhKQ0KaGVhZChDQy5Ub3RhbERhdGEpDQp0YWlsKENDLlRvdGFsRGF0YSkNCnN0cihDQy5Ub3RhbERhdGEpDQpDQy5Ub3RhbERhdGEkRmxvdy5yYXRlIDwtIGFzLmZhY3RvcihDQy5Ub3RhbERhdGEkRmxvdy5yYXRlKQ0KQ0MuVG90YWxEYXRhJENobG9yb3BoeWxsPC0gYXMuZmFjdG9yKENDLlRvdGFsRGF0YSRDaGxvcm9waHlsbCkNCg0KaW5kIDwtIHVuaXF1ZShDQy5Ub3RhbERhdGEkRF9WX1QpDQoNCmZvciAoaSBpbiAxOmxlbmd0aChpbmQpKXsNCmhpc3QoQ0MuVG90YWxEYXRhJHR1cm4uYW5nbGVbQ0MuVG90YWxEYXRhJERfVl9UPT1pbmRbaV1dLA0KICAgICAgIGJyZWFrcyA9IDUwLA0KICAgICB4bGFiID0gIlR1cm4gQW5nbGVzIiwNCiAgICAgbWFpbiA9ICIiKSANCn0NCg0KcGxvdChDQy5Ub3RhbERhdGEkRmxvdy5yYXRlLCBDQy5Ub3RhbERhdGEkdHVybi5hbmdsZSwgbWFpbiA9ICIiLCB4bGFiID0gIkZsb3cgUmF0ZSAoY20vcykiLCB5bGFiID0gIlR1cm4gYW5nbGUgKGRlZ3JlZXMpIikNCnBsb3QoQ0MuVG90YWxEYXRhJENobG9yb3BoeWxsLCBDQy5Ub3RhbERhdGEkdHVybi5hbmdsZSwgbWFpbiA9ICIiLCB4bGFiID0gIkNobG9yb3BoeWxsIChtZy9MKSIsIHlsYWIgPSAiVHVybiBhbmdsZSAoZGVncmVlcykiKQ0KDQpgYGANCg0KDQoNCmBgYHtyfQ0KZm9yIChpIGluIDE6bGVuZ3RoKGluZCkpew0KcGxvdChDQy5Ub3RhbERhdGEkdHVybi5hbmdsZVtDQy5Ub3RhbERhdGEkRF9WX1Q9PWluZFtpXV0sIGxvZzEwKENDLlRvdGFsRGF0YSR2W0NDLlRvdGFsRGF0YSREX1ZfVD09aW5kW2ldXSksDQogICAgICAgICAgICB4bGFiID0gIlR1cm4gQW5nbGVzIiwNCiAgICAgeWxhYiA9ICJWZWxvY2l0eSIsDQogICAgICAgICAgbWFpbiA9ICIiKSANCn0NCg0KYGBgDQoNCg0KYGBge3J9DQpzYXZlLmltYWdlKCJ+L1Bvc3QtZG9jL0RhdGEvVG90YWwgTWVyZ2VkIERhdGEgRmlsZS5SRGF0YSIpDQoNCmBgYA0KDQoNCg0K